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Abstract. The High Mass X-ray Binary Pulsar SAX J2103. 5+4545 has been observed with INTEGRAL several 
^ ' times during the last outburst in 2002-2004. We report a comprehensive study of all INTEGRAL observations, 

allowing a study of the pulse period evolution during the recent outburst. We measured a very rapid spin-up 
episode, lasting ~130 days, which decreased the pulse period by ~1.8 s. The spin-up rate, P~— 1.5xl0~ 7 s s" 1 , is 
the largest ever measured for SAX J2103. 5+4545, and it is among the fastest for an accreting pulsar. The pulse 

in ■ 

profile shows evidence for temporal variability, apparently not related to the source flux or to the orbital phase. 
The X-ray spectrum is hard and there is significant emission up to 150 keV. A new derivation of the orbital 
period, based on RXTE data, is also reported. 
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, 1. Introduction and spin periods in Be/X-ray binaries (which would imply 

an orbital period ^190 days). 

<3 • SAX J2103. 5+4545 is a transient Be High Mass X-ray „, , ... , , . , 

* • /TTw^ns , /! -, <^ont -Lhe absorption m the direction oi the optical coun- 

>• , Binary (HMXRB) pulsar (pulse period ^358s) discov- . ,. , ,. . r „ . , . nrtnA ^ 

, . , , f ^. , , „ ,„ mn , 7 , , terpart indicates a distance of 6.5 kpc (Reig et al., 2004), 

k ^ ■ ered with the Wide Field Cameras (WFC) on-board , ,,. - „ „, m , \ , / 

r>^ _ „ . „ . . _ / ' . doubling the previous estimate ol 2 kpc (Baykal et al., 

r> , BeppoSAX during an outburst m l H ebruary-beptember „„„ nN , , , , c . , ^ ^ 

H ^^/tt i, , i ,nn„\ a i • , •; ^ i 2002) based on a model lor the spin-up produced by ac- 

_a : 1997(HuUeman et al. 1998 A peak intensity of 20 mCrab ^ g & » 1 

was reached on 1997 April 11. Ihe X-ray spectrum in the 

2-25 keV energy band was fit by a power-law with a pho- A transition from the spin-up phase to the spin-down 

ton index r=1.27±0.14 modified by photo electric absorp- re S ime was observed with RXTE during the decline in 

tion at lower energies (JV H =3.1 x 10 22 cm" 2 ). the X " ra y flux of the 1999 outburst (Baykal et al., 2002). 



An orbital period of 12.68 days and an orbital eccen- 
tricity, e, of 0.4+0.2 have been measured with the Rossi 



During the initial part of the outburst the source exhib- 
ited a spin- up phase (the pulse period decreased by ^0.9 s 

. _ , ,,,, , , , , in 150 days), but when the flux dropped the pulse fre- 

X-ray liming Explorer (RX IB) during a second outburst , , „ , , , ,. 

f , , , , , ,. ° /T , , , , quency saturated and, as the flux continued to decline, a 
lasting about 14 months observed m 1999 Baykal et ah, , . , , , . , . . , 

, weak spm-down phase started, this correlation between 

spin-up rate and X-ray flux indicates the formation of an 

The likely optical counterpart is a BOVe star (V=14.2; accretion disk during the periastron passage. The X-ray 

Reig & Mavromatakis, 2003; Reig et al. 2004), making spectrum, described by a power-law with a cut-off around 

SAX J2103. 5+4545 the Be/X-ray binary with the shortest g keV, displayed an iron emission line at 6.4keV (Baykal 

orbital period known. SAX J2103. 5+4545 does not follow e ^ a j_, 2002). 

the correlation found by Corbet (1986) between the orbital „ . V n^n^^ T i ^ r i • 

v ' During a XMM-JNewton observation performed m 

January 2003, Inam et al. (2004), detected for the first 

Send offprint requests to: L. Sidoli (sidoli@mi.iasf.cnr.it) time a soft blackbody component with a temperature of 
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1.9 keV and radius 0.3km (assuming a 3.2 kpc distance), 
probably produced in the polar cap of the neutron star, 
as well as a transient quasi-periodic oscillation at 22.7 s. 
Besides the blackbody component, the 1-20 keV spectrum 
(using XMM-Newton and RXTE data) was fit with a 
power-law with a cut-off at 7.89 keV and an e-folding 
energy of 27.1 keV, together with an emission line at 
6.42 keV. The equivalent width of this cold iron line was 
found to depend on the spin phase. Using RXTE data 
(MJD 52611.48 - MJD 52668.90), Inam et al. (2004) de- 
rived a pulse period of 354. 7940±0. 0008s (at the epoch 
of the XMM-Newton observation, MJD 52645.85) and a 
57.5 days average spin-up rate of (7.4+0.9) xlO -13 Hz s _1 . 

Analyses of part of the INTEGRAL observations of 
the source have been reported by Lutovinov et al. (2003), 
Sidoli et al. (2004), Blay et al. (2004), Filippova et al. 
(2004) and Falanga et al. (2005). In this paper we report 
on the timing and spectral analysis of all INTEGRAL ob- 
servations performed between December 2002 and July 
2004. 

2. Observations 

The ESA INTEGRAL observatory, launched in October 
2002, carries 3 co-aligned coded mask telescopes: the im- 
ager IBIS (Ubertini et al. 2003), which allows high-angular 
resolution imaging in the energy range 15keV-10MeV, 
the spectrometer SPI (Vedrenne et al. 2003; 20 keV- 
8MeV) and the X-ray monitor JEM-X (Lund et al. 2003; 
3-35 keV). IBIS is composed of a low-energy CdTe detec- 
tor ISGRI (Lebrun et al. 2003), sensitive in the energy 
range from 15keV to lMeV, and a Csl detector (PICsIT; 
Labanti et al. 2003), designed for optimal performance 
at 511 keV, and sensitive in the 175keV-10MeV energy 
range. 

The INTEGRAL satellite performs regular scans of the 
Galactic Plane (GPS) every 12 days as part of the Core 
Programme to monitor the timing and spectral properties 
of known X-ray sources, discover new transient sources 
and map the diffuse emission from the Galactic plane 
(Winkler et al. 2003). 

From December 2002 to July 2004, SAX J2103. 5+4545 
has been observed several times. Among these observa- 
tions (170 pointings), some have been performed dur- 
ing the satellite performance and verification phase 
(December 2002), during the GPS (January 2003-July 
2004), or during guest observer observations now in the 
public archive (May 2003). We have performed a com- 
prehensive study of all these observations (see Table , 
concentrating in particular on the timing properties of 
SAX J2103. 5+4545. 

The data reduction has been performed using the 
OSA 4.2 release of the INTEGRAL analysis software, 
with the corresponding response matrices. Each point- 
ing has a typical exposure of ~2-3ks. Only observations 
where the SAX J2103. 5+4545 position offset was less than 
9° have been considered. Only for a few of these point- 
ings also JEM-X data were available, due to the smaller 



field of view of this instrument (offset less than 3.5°). 
SAX J2103. 5+4545 has been detected in the energy range 
20-40 keV with IBIS/ISGRI in 113 pointings (see Tabled). 

The 2-10 keV lightcurve of SAX J2103. 5+4545 ob- 
tained with the RXTE All Sky Monitor (ASM) from 1997 
to 2004 is shown in Fig. d The times of the two previous 
outbursts observed with BeppoSAX (in 1997) and RXTE 
(in 1999), as well as the epoch of the INTEGRAL obser- 
vations reported here, are indicated. 



Table 1. Summary of the INTEGRAL observations of 
SAX J2103. 5+4545. The observations have been grouped 
together for brevity. The 1st and 2nd columns list the 
Start and Stop Time of each group of observations, col- 
umn n.3 reports the number of ISGRI pointings where 
the source is within 9° of the pointing direction, column 
n.4 reports the number of ISGRI/IBIS pointings where 
SAX J2103. 5+4545 has been detected in the energy range 
20-40 keV, while the last column shows if a search for pe- 
riodicity has been succcsfully carried out. 



Start Time 


End Time 


Number 


Detections 


Timing 


(MJD) 


(MJD) 


of obs. 






52617.86 


52619.28 


4 


3 


N 


52620.74 


52620.93 


2 





N 


52629.63 


52631.92 


58 


56 


Y 


52636.66 


52637.87 


21 





N 


52638.46 


52640.07 


2 





N 


52653.32 


52653.40 


2 


1 


N 


52722.85 


52722.93 


3 





Y 


52737.03 


52737.12 


3 


1 


Y 


52746.10 


52746.18 


3 


3 


Y 


52761.26 


52762.42 


38 


38 


Y 


52770.05 


52770.10 


2 


2 


Y 


52782.06 


52782.14 


3 


3 


Y 


52835.9 


52836.03 


3 


1 


N 


52985.41 


52985.49 


3 





N 


52994.38 


52994.47 


3 





N 


53019.37 


53019.45 


3 





N 


53045.61 


53045.69 


3 





N 


53090.96 


53091.00 


1 





N 


53102.48 


53102.56 


3 





N 


53114.16 


53114.24 


3 





N 


53126.09 


53126.14 


2 


2 


Y 


53137.88 


53137.93 


2 





N 


53164.97 


53165.05 


3 





N 


53188.96 


53189.03 


2 





N 
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Fig. 1. ASM/RXTE lightcurve (8-days average) of 
SAX J2103. 5+4545 from the outburst in 1997 to the end 
of August 2004. Three distinct outbursts are evident, and 
have been marked as "SAX 1997" (time of the source 
discovery with BeppoSAX, Hulleman et al. 1998), "XTE 
1999" (observations with XTE reported by Baykal et al. 
2000), "2002-2004" outburst observed with INTEGRAL. 



3. Results 

3.1. Timing Analysis 

In order to perform the timing analysis, we selected events 
with a Pixel Illumination Function (PIF, pixel fraction 
illuminated by the source) larger than 0.5. IBIS/ISGRI 
events have been collected in the energy band 20-40 keV, 
where the best statistics is available, while for JEM-X the 
range 3-30 keV has been used. 

Arrival times have been corrected to the Solar System 
barycenter and for the Doppler delay due to the mo- 
tion in the binary system. To perform this orbital 
correction, we adopted the ephemeris zero point from 
Baykal et al. (2000), but used an updated orbital pe- 
riod, P or b=12. 670+0. 005days, that we calculated with 
epoch folding of all the available RXTE/ASM data (see 
Fig. Eland Fig- EJ ■ F° r °ur analysis, based on ASM dwell- 
by-dwell data, we used the epoch folding method de- 
scribed in Larsson (1996). The error estimate is based 
on Monte Carlo simulations and on the analysis of pul- 
sations injected into the ASM lightcurve. From the Monte 
Carlo simulations we estimate the statistical error to be 
±0.003 days. The slightly higher error estimate above in- 
cludes also effect of the irregular source variability on the 
period determination. This latter effect is estimated from 
the injected pulsation. 

Since the short duration of each single INTEGRAL 
pointing (~2-3ks) leads to a large uncertainty in the de- 
termination of the pulse period (typically about 10 s, even 
in presence of good statistics), we grouped the observa- 
tions where a good estimate of the pulse period could be 




12.40 



12.50 12.60 12.70 

FOLDING PERIOD (Days) 



12.80 



12.90 



Fig. 2. x 2 distribution using all avilable RXTE-ASM 
data. It peaks at 12.670+0.005 days, our improved deter- 
mination of the orbital period for SAX J2103. 5+4545. The 
solid line marks the data and the dashed line the fitted 
model (see text). 
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Fig. 3. SAX J2103.5+4545 RXTE-ASM lightcurve folded 
at the updated orbital period of 12.670+0.005 days. 



achieved in 7 separate data-sets consisting of the pointings 
performed during the same satellite revolution. 

In Table [5] we list the start and end time of each 
group of observations, together with the source average 
IBIS/ISGRI count rate in the 20-40 keV range. In data- 
sets n.l, n.4 and n.5 also JEM-X data are available, and 
have been used in the timing analysis. 

We searched for the SAX J2103. 5+4545 spin period by 
using epoch folding techniques: we produced x 2 distribu- 
tions versus trial period, and then we fit them in order 
to derive the best-fit period for each data-set, and its as- 
sociated uncertainty (Leahy 1987). An example of a \ 2 
distribution is shown in Fig. 

The derived pulse periods are reported in Tabled Our 
data-sets n.l and n.4 include observations partly reported 
by Blay et al. (2004), Filippova et al. (2004) and Falanga 
et al. (2004). 
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Table 2. Results of the SAX J2103. 5+4545 timing analysis. The 2nd and 3rd columns list the Start and Stop Time 
of each data-set, the 4th the averaged IBIS/ISGRI (20-40 keV) source count rate in these epochs, the 5th and 6th 
columns list the pulse periods measured in the seven INTEGRAL data-sets. In data-sets n.2, 3, 6 and 7 the source 
was outside the JEM-X field of view. In data-set n.l only the pulse period estimated with JEM-X is reported due to 
several telemetry gaps present in ISGRI lightcurve during this satellite revolution, which hampered a proper timing 
analysis 



Data-set 
ID. 


Start Time 
(MJD) 


End Time 
(MJD) 


IBIS/ISGRI 

(s- 1 ) 


Pulse Period (s) 
IBIS/ISGRI JEM-X 


1 


52629.91 


52631.84 


5.16+0.20 




355.09 ±0.03 


2 


52722.85 


52722.90 


4.28+0.18 


356. ±3 




3 


52746.10 


52746.18 


6.07+0.25 


355.5 ±1.0 




4 


52761.26 


52762.42 


4.07+0.08 


353.32 ±0.03 


353.47 ±0.15 


5 


52770.05 


52770.10 


8.62+0.27 


351.8 ±1.3 


355.4 ±2.3 


6 


52782.06 


52782.12 


6.45+0.28 


351.3 ±2.3 




7 


53126.09 


53126.14 


4.34+0.28 


355. ±3 






345 350 355 360 

PERIOD (s) 



Fig. 4. Distribution of \ 2 versus trial period for data-set 
n.4 (see Table |3J) in the energy range 20-40 keV 



A period difference of ~1.77s in ~130.5days is evi- 
dent from the comparison of results from data-sets n.l 
and n.4, which have the smallest uncertainty (Table 
This spin-up is significantly higher than that measured 
during previous outbursts of this source. 

The pulse period evolution during the INTEGRAL 
observations is shown in Fig. [7J where also the period 
measured by Inam et al. (2004) with RXTE has been 
shown for completeness. A linear fit to all these data in- 
dicates an average spin-up with a rate of P=— (1.487 ± 
0.030) xl0" 7 ss- 1 (x 2 =21.3, 8 degrees of freedom, dof). 

A few examples of the folded light curves (background 
subtracted) obtained in different observations with ISGRI 
and JEM-X are shown in Fig. and Fig.El respectively. 



3.2. Spectral Analysis 

We used data from all INTEGRAL observations to derive 
a spectrum with a very high signal to noise ratio. Since 
the source exhibits no significant spectral evolution during 
the outburst, this can be safely done. 

The average source spectrum from 5 to 200 keV, has 
been investigated by simultaneously fitting JEM-X (5- 
20keV) together with IBIS/ISGRI (20-200 keV), and SPI 
(20-200 keV) spectra, achieving net exposure times of 
40 ks, 290 ks and 130 ks, respectively for the three in- 
struments. An additional uncertainty of 5% was added 
quadratically to the statistical errors to account for sys- 
tematic uncertainties in instrumental responses. Free rel- 
ative normalizations between the three instruments were 
included. Since the energy range (5-200 keV) is not sen- 
sitive to the interstellar column density, we fixed it at 
3xl0 22 cm~ 2 (Hulleman et al., 1998). 

A single power-law is not able to describe the broad- 
band spectrum (reduced ~x 2 >2). The best-fit (x 2 = 121.9 
for 129 d.o.f.) is an absorbed power-law with a high 
energy cut-off (see Fig. [SJl, with the following parame- 
ters: photon index T=1.53 l"o!o7> cu t-°ff energy E c =19 
t^keV and e-folding energy E f =32±3keV. The 5-200 
keV observed flux is ~1.05xl0 _9 ergcm -2 s _1 (based on 
the JEM-X response matrix). The source flux in the 
2-10 keV energy range, corrected for the absorption, is 
^4.2xl0~ 10 ergcm~ 2 s _1 , which translates into a lumi- 
nosity Lx=2xl0 36 ergs -1 (2-10 keV, for a distance of 
6.5 kpc). 

4. Discussion 

We have studied all the available INTEGRAL observa- 
tions of SAX J2103. 5+4545, obtained from December 2002 
to July 2004. They cover part of the latest source out- 
burst, which is the longest and brightest observed from 
this source since its discovery in 1997. 
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Fig. 5. Few examples of SAX J2103. 5+4545 pulse profiles from INTEGRAL observations where a different shape is 
evident. From top to bottom, from the left to the right, profiles are from data-sets n.2, n.5, n.6 and n.7, all in the 
energy range 15-40 keV. The zero phase is arbitrary. The pulsed fraction, defined as the ratio of maximum minus 
minimum to maximum, is ^33% for data-set n.2, while ^47%-50% for the other three lightcurves 



SAX J2103.b+4b4b 
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Fig. 6. Pulse profiles at soft energies (3-30 keV) with JEM-X, from data-set n.l (left panel) and n.4 (right panel). 



To perform a proper timing analysis of the 
INTEGRAL observations, we re-determined the orbital 
period of the binary system using ASM/XTE public data, 
obtaining a refined P or b of 12.670±0.005days (to be com- 
pared with P or b of 12.68±0.25days, measured with XTE 
during the 1999 outburst by Baykal et al., 2000). 

The SAX J2103. 5+4545 pulse profiles observed with 
INTEGRAL show evidence for temporal variability over 



the same energy range (15-40 keV), where a single asym- 
metric pulse profile seems to evolve to a double-peaked 
profile in some of the observations. On the other hand, 
in the lightcurves extracted at softer energies (3-30 keV) 
there is no clear evidence for a secondary peak. The 
INTEGRAL observations unfortunately do not cover the 
entire orbital phase; in particular, the data-sets n.3, n.5, 
n.6 and n.7 sample roughly the same orbital phase (these 
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2 — INTEGRAL Obs. 

= III I | 

2 1.5r n 




MJD 

Fig. 7. Upper panel: A close-up view of the ASM/RXTE lightcurve (8-days average) during the INTEGRAL 
observations. The epochs of the INTEGRAL observations have been marked with vertical lines Lower panel: 
SAX J2103. 5+4545 pulse period evolution (triangles mark JEM-X results; squares IBIS/ISGRI results). The dia- 
monds indicate the period measured with XTE by Inam et al. (2004). A linear fit to all the data is reported (solid 
line) together with the 1-a uncertainty (dashed lines). 



are indeed GPS observations, which are performed ev- 
ery 12 days, a value close to the source orbital period). 
Thus, apparently the complex variability of the pulse pro- 
file does not depend on the orbital phase. Also an obvious 
dependence on the source luminosity seems to be unlikely, 
since we observe different pulse periods from observations 
at similar fluxes (e.g., data-sets n.2 and n.7). The com- 
plex variability of the pulse profiles is a common feature 
among the class of accreting binary pulsars, also over the 
same energy range (White, Swank & Holt, 1983; Bildsten 
et al., 1997). In some sources the pulse shape variability 
shows a clear dependence on the source luminosity (first 
observed in EXO 2030+375, Parmar et al., 1989) or on 
the energy range (e.g. Vela X-l, Kretschmar et al, 1997). 
A simple geometric model with two pencil beams coming 
from the magnetic poles is too simple to explain the pulse 
profiles in most pulsars, and a more complex emission 
pattern or the inclusion of scattering effects in the trans- 
fer through highly magnetized plasmas are needed (e.g., 
Nagel, 1981). A computation of pulse profiles of medium- 
luminosity (Lx ^10 36 erg s _1 ) binary X-ray pulsars has 
been performed by Kraus et al. (2003), who found that 



energy-dependent peaks in the profiles are mainly due to 
energy-dependent relative importance of the halo (which 
forms around the accretion funnel where the neutron star 
surface is irradiated) and the column contributions to the 
observed flux. It is also possible that the temporal vari- 
ability of the pulse profile we observe in the same energy 
range can be due to time-dependent emission pattern, or 
to changes in the opacity of the magnetized plasma where 
the radiation propagates. 

The pulse period evolution measured with INTEGRAL 
shows clear evidence for a spin-up episode, lasting at 
least ~130days from MJD 52629.9 (data-set n.l) to 
MJD 52761.3 (data-set n.4), where the pulse period de- 
creased by ~ 1.77 s. 

Considering the overall pulse period history during the 
latest outburst (Fig.EJ, a linear fit results in a spin-up rate 
P=-(1.487 ± 0.030) x 10~ 7 s s" 1 . 

We cannot exclude a more complex evolution of the 
pulse period than a simple steady spin-up over more than 
1.3 years. In particular, it is likely that the spin-up trend 
induced by the accretion does not extend until the epoch 
of the data-set n.7. 
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Fig. 8. Left panel: SAX J2103. 5+4545 counts spectrum (5-200 keV) extracted from JEM-X (5-20 keV), IBIS/ISGRI 
(20-200 keV, upper counts spectrum) and SPI (20-200 keV, lower spectrum) instruments. The best-fit model (his- 
togram) is an absorbed power-law with a cut-off at high energies (see text for the best-fit parameters). The residuals 
displayed in the lower panel are in units of standard deviations. Right panel: SAX J2I03. 5+4545 unfolded spectrum 
together with the best-fit model (solid line). 



A variation in the accretion rate can modify the 
pulse period derivative. A correlation of the spin-up 
rate with the X-ray flux has been indeed found for 
SAX J2I03. 5+4545 by Baykal et al. (2002) during the 1999 
outburst, as predicted by the Ghosh and Lamb (1979) 
model for the torques exerted on the neutron star during 
the accretion through a disk. The existence of this cor- 
relation is indeed indicative of the formation of an accre- 
tion disk around SAX J2103. 5+4545 during that outburst. 
The last pulse period measured with INTEGRAL (data- 
set n.7) has been determined after a long faint state last- 
ing about 300 days (Fig. upper panel). The reduced ac- 
cretion rate could have increased the SAX J2I03. 5+4545 
pulse period (a spin-down phase followed the spin-up trend 
during the 1999 SAX J2103. 5+4545 outburst, Baykal et 
al. 2002), although the large uncertainty in the pulse pe- 
riod at that epoch (P=355±3s) does not permit to draw 
firm conclusions. 

The formation of an accretion disk in 
SAX J2103.5+4545 during the outburst in 1999 has 
been established from RXTE data (Baykal et al. 2002). 
Furthermore, the fact that our estimate of the spin-up 
rate (the highest ever measured in this source) has 
been measured during the brightest outburst observed 
from SAX J2103. 5+4545, indicates that the correlation 
between luminosity and spin-up is still valid. It is thus 
tempting to discuss our results on the source spin-up rate 
in the framework of the Ghosh and Lamb (1979) theory. 

In Fig. we compare the measured pulse period 
derivative of ~— 1.49xl0 _7 ss _1 with the theoretical re- 
lation between the spin- up rate of a 1.4M neutron 
star and the quantity PL 3 / 7 , where P is the pulse pe- 
riod in seconds, and L is the source luminosity in units 
of 10 37 ergs _1 (Ghosh & Lamb, 1979). Thus, assuming 
that the distance of 6.5 kpc is correct, a magnetic mo- 
ment of ~1.6x 10 30 Gauss cm 3 is needed to explain the 
spin-up rate. This translates into a magnetic field of 



~10 12 Gauss, which would imply a cyclotron line at ener- 
gies around ^10-20 keV, not observed in the source spec- 
trum. The absence of lines could suggest a much higher 
neutron star magnetic field, ~10 13 Gauss, but this would 
imply a smaller distance, around 5 kpc. A distance of 
5 kpc also agrees with the maximum spin-up trend mea- 
sured with XTE during the 1999 outburst (Baykal et 
al., 2002): a pulse period derivative of ~8.9x 10 -8 s s _1 
(and pulse period P=358.166s) for a source flux of 
^6x 10~ 10 erg cm -2 s _1 can be explained with a neutron 
star of 1.4 M Q and a magnetic moment ~10 31 Gauss cm 3 , 
located at a distance of about 5 kpc. 

However, note that there are uncertainties in the di- 
mensionless angular momentum adopted in the accretion 
torque theory (Li & Wang, 1996), which could significantly 
modify the estimate of the source distance derived with 
this method (see, e.g., Tsygankov & Lutovinov (2005), on 
the binary pulsar KS 1947+300). Another uncertainty is 
due to the source flux we adopted here, which is not the 
bolometric flux and has been measured from the average 
INTEGRAL spectrum in the energy range 5-200 keV. 

The spin-up rate measured with INTEGRAL from 
SAX J2103. 5+4545 is the largest experienced by this 
source (likely because of the higher brightness of the latest 
outburst) and is among the largest found for an X-ray bi- 
nary pulsar. The frequent monitoring (as that performed 
with INTEGRAL during the GPS) of accreting pulsars is 
essential in the determination of the pulse period evolu- 
tion, and crucial to test accretion theories. 
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